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Abstract 

We study the structures of hybrid stars with leptons at finite temperature under beta equiUbrium. 
For the quark phase, we use the three flavor Nambu-Jona-Lasinio (NJL) model. For the hadron 
phase, we adopt nuclear equation of state (EOS) by Shen et al.. This EOS is in the framework 
of the relativistic mean field theory including the tree body effects. For the hadron-quark phase 
transition, we impose the bulk Gibbs construction or the Maxwell construction to take into account 
uncertainties by finite size effects. We find that the pure quark phase does not appear in stable 
star cores in all cases. With the phase transition, the maximum masses increase ~ 10% for high 
lepton fraction. On the contrary, without the transition, they decrease ~ 10%. We also find that, 
in the NJL model, the lepton fraction is more important for structures of unstable stars than the 
temperature. This result is important for many astrophysical phenomena such as the core collapse 
of massive stars. 
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I. INTRODUCTION 



The possible existence of quark matter in compact stars has been proposed l|,l2j. It should 
be noted that the uncertainty of EOS is always a big problem in the study of compact star 
structures. The recent observations imply the existence of large mass compact objects 4, 

The effects of quark matter on various astrophysical phenomena have been studied ex- 
tensively. For example, cooling effects on compact star evolutions have been studied in 



Ref. 
radiation 



9, 



11 



10 [. Other examples include the effects of quark matter on gravitational wave 
12, 13], neutrino emissions fl^. [l5|. rotational frequencies 16], and the max- 
imum energy release during the transition from neutron stars to quark stars |l7, [l^, etc.. 
However, most of them were based on the MIT bag model. 

On the other side, many researchers have proposed EOS's based on finite temperature field 
theory in beta equilibrium. Among these, the Nambu-Jona-Lasinio (NJL) model is widely 
known. The difference between the NJL model and the MIT bag model on the structure 
of compact stars has been pointed out in Ref. 2^. Furthermore, the NJL model has 
two notable points. One of them is that the NJL model implies a color-superconducting 
state, which effects cooling evolution on compact stars with quark matter 2l|. The other is 
that it can show chiral phase transition. We can calculate the spontaneous chiral symmetry 
breaking by the NJL model. However, baryons in the framework of the NJL model do not 
reflect experimental result well. This is a big problem in the adaptation of the NJL model 
for realistic astrophysical phenomena. 

To study realistic compact stars including exotic matter, Burgio and Plumari (2008) have 
set two kinds of EOS's for hadron and quark phases [l^. Hadron EOS in their study was 
based on Bruekner-Bethe-Goldstone many body theory, which was considered appropriate 



for neutron matter 



22l ]. Furthermore, they adopted the NJL model for the quark phase. 



and imposed the Maxwell construction on the hadron-quark phase transition. However, 
there are many uncertainties for the hadron-quark phase transition. Assuming the quark 
deconfinement transition to be of first order, it causes a thermodynamical instability and 
the mixed phase appears around the critical density. The properties of the mixed phase 
depend on electromagnetic and surface contributions to the energy. These contributions are 
sometimes called "finite-size effects" . The quantitative analysis of these finite size effects has 
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provided following result: EOS's for the phase transition become similar to the ones under 



the bulk Gibbs construction for weak sur: 
construction for strong surface tension 23 



ace tension, and to the ones under the Maxwell 



24 



25|. 



In the present paper, we adopt a widely accepted EOS in astrophysics for the hadron 
phase. The EOS is based on the relativistic mean field (RMF) theory including many body 



effects 



261 ] ■ For the quark phase, we also adopt the NJL model. We impose the bulk Gibbs 



or the Maxwell construction on the phase transition to take into account uncertainties of 
finite size effects. For realistic proto-neutron star (PNS) structures, the finite temperature 
effects and neutrino trapping are important. For this reason, in this paper, we also consider 
the finite temperature case with neutrinos as well. 

This paper is organized as follows. In Sec. II, we outline our EOS's. In Sec. Ill, we present 
numerical result for PNS structures. Sec. IV is devoted to the conclusion and discussion 
where we discuss the astrophysical meaning of our result. 



II. EQUATION OF STATE 

A. Equation of state for quark phase — Three flavor Nambu-Jona-Lasinio model 



In a moderately dense system, for example inside compact stars, the quark matter may 
exist with s— quarks. 

The Lagrangian density of the NJL model in the three flavor system is written as 



where 



C = q{i ^ - ■m)q + C4 + Cq, 



fc=0 



-K det(gi(l + 75)5^) + det(gi(l - 75)9^^ 



(1) 
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Here, i and j denote the flavor indices, whereas q and denote quark field and current mass 
matrix, respectively. The Lagrangian density components, £4 and Cq respectively generate 
four- leg and six- leg interaction. The six- leg interactions, Cq, come from anomaly. In 

the chiral limit, mj = , the Lagrangian density has SU{3)l x SU{3)fi x SU{3)c x f/(l)v 
symmetry. 

With the mean field approximations, the Lagrangian density is given as 

C-MFA = ^ qf{i^-{mf + ^f^^))qf-U, (4) 

where 

Si,s = -i4:Gs{qiqi) -2K{qjqj){qkqk)), (5) 
U= Y.{2GMqif)-^K{q^qu){md){qsqs), (6) 



l=u,d,: 



for f ^ j ^ k. 

Therefore, the thermodynamic potential is 



O T 

- = --\nZ (7) 
V V ^ ^ 



f=u,d,s ' 
1 , 1 



+ ' ln(l + e-'^^/^P)) + ^ ln(l + e'^^^P); 
+U, (8) 



where ^/(p) = y^p^ + M|, = Ef{p) ± /// , and M/ = m/ + S/,s- 

Since the NJL model is nonrenormalizable, it is necessary to introduce a cutoff in the 
momentum integration. In this study, we use the three-dimensional momentum cutoff and 
we use the parameter set that is obtained in Ref. The cutoff A is 0.6023 GeV. The 

coupling constant of four quark interaction is defined as GgA^ = 1.835, while that of six 
quark interaction is defined as KA^ = 12.36. The current masses of u— and (i— quarks are 
fixed at 5.5 MeV and that of s— quark is fixed at 140.7 MeV. 

The quark number density Ug is given as Ug = + Ud + ris = Sub with nj = {q^q/} 
where is the baryon number density. Moreover, the chiral condensate (q/qj) satisfies the 
stationary condition dQ/d{qfqf) = 0. 
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FIG. 1: The left panel shows the pressure versus the baryon number density ns under beta 
equilibrium and charge neutrality in the NJL model. The right panel shows the effective quark 
masses versus hb- Here, we normalized by the nuclear saturation number density no; no = 
0.17/?7i~^. We do not show the effective masses of d— quark in this figure, because they are almost 



the same values as the ones of ?i— quark 



2C|. 



Inside of compact stars, we must impose beta equilibrium and charge neutrality. There- 
fore, we rewrite the chemical potential of u—, d— and s— quarks as 



2 A'-l'e ) 



(9) 
(10) 



where /i^ is the effective baryon chemical potential corresponding to the neutron chemical 
potential. We assume that neutrinos are trapped because it is reasonable for PNS's in the 
process of core collapse of massive stars. 

Our equation of state in the NJL model is shown in FIG. [T] We set the temperature as 
T = — 30 MeV in this figure, because this is the expected temperature range in compact 
stars 281]. For the total lepton fraction, we adopt Yi = 0.1 — 0.4 for the same reason. 
The lepton fraction is more important in deciding the stiffness of EOS's than temperature 
because of the following reason: A high lepton fraction provides a high electron fraction 
which suppresses the number of other negatively charged particles, such as s— quarks, due 
to the electric charge neutrality condition. This implies that at high lepton fraction case, 
being the strangeness of content small, the chiral restoration for s— quark is suppressed as 
shown in the right panel of FIG. [H Hence, EOS's for high lepton fraction become stiff 20] . 



However, the above discussion does not apply to the MIT bag model, because the MIT bag 
model does not contain the effects of chiral restoration. The details are shown in Menezes 
et al.(2005) 20| or Burgio et al. (2003, 2008) 

An important thing for our EOS is the mass of the s— quark and the number density 
which primarily depend on the chiral restoration, the electric charge neutrality, and the 



Gibbs/Maxwe 



our result 
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1 construction. However, effects of the color superconductivity may change 



30, 



31 



32 



33, 



34 



35l | . Studies of such effects are out of scope of this work. 



In the NJL models, the bag pressures are not free parametric constants as in the MIT 
bag model. 



B. Equation of state for Hadron phase and Mixed phase 

For the EOS of the hadron phase, we adopt the nuclear EOS based on the RMF theory 



developed by Shen et al. (Shen EOS) 26|]. This EOS includes three-body effects and has been 



data of masses and radii of stable and unstable 



261]). The range of baryon density is from 10^'^ g 



constructed to reproduce the experimenta 
nuclei (see references in Shen et al.(1998) 
cm~^ to 10^^'^ g cm~^. 

Pressure as a function of baryon number density is shown in FIG |2l Contrary to the NJL 
model, this EOS becomes slightly soft for high lepton fraction: e.g. the pressure at Yi = 0.4 
and T = 30 MeV is about 10 % lower than the one at = 0.1 and the same temperature. 
The explanation for this is that, under charge neutrality condition, a high lepton fraction 
decreases the neutron number density which is the main component of a repulsive nuclear 
interaction above the saturation density. 

For the phase transition, there are many uncertainties which come from finite size effects 
as discussed in Sec. I. In the past studies, it has been shown that EOS's become similar to 



the ones under the Maxwell construction for strong sur 



23|, 



'ace tension, and to the ones under 



24 



251]. Hence, to take into account 



the bulk Gibbs construction for weak surface tension 
such uncertainties, we calculate EOS's under the bulk Gibbs construction and the Maxwell 
construction, respectively. 

First, we describe the bulk Gibbs construction. In this construction, we consider chemical 
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FIG. 2: Same as the left panel of FIG. [H but for the nucleoli EOS by Shen et al. 



equilibrium at the hadron-quark interface as well as in each phase: 



26| 



Hp + He 

Hn 



Hn ~\~ Hl^e 1 

flu + 



(11) 
(12) 
(13) 
(14) 



Here, fip fin are the proton chemical potential and the neutron chemical potential, 
respectively. We assume that fin is the same as /i^ in equation ([9]) and ( fTOl) . 

The bulk Gibbs constructions impose the following other conditions for phase equilibrium: 



pQ ^ pH 



(15) 



Here, and are the temperatures in the quark phase and in the hadron phase, respec- 
tively. Similarly, P*^ and denote the pressures in those phases. 

Besides, we consider baryon number conservation law in the mixed phase as. 



flB = X'lT'B,, 



x)nB,H 



where nB,nB,Q,nB,H, and x are respectively the total baryon number density, the baryon 
number density in hadron matter, the effective baryon number density in quark matter, and 
the volume fraction which shows quark matter volume divided by the total volume. Here, 
nB,Q is one third of the total quark number density. 
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Furthermore, we adopt the electric charge neutrahty as follows; 



Yens = X yp,Q nB,Q + (1 - x) ^p,/f ^b^h ■ 

Here Ye is the electron fraction, whereas Yp^g and Yp^n denote the fraction of positive particles 
per baryon in quark matter, and in hadron matter. 

Schertler et al. 36| have calculated the mixed phase under the bulk Gibbs constructions. 
Their EOS for hadron phase was also based on the RMF theory, and their EOS for quark 
phase was also the NJL model. However, their calculations were only for zero temperature 
and did not include neutrinos. 

Under the Maxwell construction, lepton chemical potentials are not considered in the 
phase equilibrium. The condition of chemical equilibrium for the phase transition is only 
equation f[T^ . The Maxwell construction imposes equation f[T5]) . too. In the following, we 
replace T^(T«) and P^(P«) by T = = T« and P = = pQ. 

Pressure as a function of baryon density is shown in FIG. [3l In this figure, we set the 
temperatures and the total electron fractions as T = — 30 MeV and Yi = 0.1 — 0.4 as 
in subsection A. At these temperatures and lepton fractions, hybrid stars are stable with 
low central densities, and have maximum masses as we show later in figures of mass-radius 
relations. Under the Gibbs construction, the density regions of mixed phase are wide; e.g. 
from 1.19 riQ to 8.40 uq for T = MeV and Yj = 0.1, where riQ is the nuclear saturation 
number densit y g iven by 0.17 fm~'^. Such a wide range of mixed phase was also found by 



Schertler et al [36|]. Under the Maxwell construction, the EOS's have density jumps from the 
purely hadronic phase to the pure quark phase; e.g. from 2.37 uq to 3.54 uq for T = MeV 
and Yi = 0.1. We also find that the quark phase appears in low density with increasing 



temperature, as shown in Burgio et al 



19|. 



III. STRUCTURE OF HYBRID STARS 



We calculate the spherical structures of hybrid stars by Tolman-Oppenheimer-Volkoff 
(TOV) equation using the EOS's given in Sec. II. These PNS's have a temperature range of 
T = — 30 MeV, and a lepton fraction range of Y; = 0.1 — 0.4. 

In FIG. mini we display mass versus central density and mass versus radius for compact 
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FIG. 3: Pressure as a function of the baryon density normalized with the nuclear density under 
the bulk Gibbs construction (left panel) and the Maxwell construction (right panel). 



stars. We assume the dynamical stability condition as 

^ > 0, (16) 

where M and ub^c are the stellar mass and the baryon number density of stellar core, 
respectively. 

In this paper, we calculate the structures of PNS's at constant temperature for simplic- 
ity, not as isentropic matter. The former situation is relevant to the neutrino transparent 
case (small Yj) whereas the latter, to the supernova stage (large Yi) |37l |. 



A. Stable stars 

In this subsection, we discuss the stars satisfying condition f|T6|) for each EOS. We find 
that stable stars do not have pure quark phase in their cores for all cases. 

First, we discuss the stable stars with the hadron-quark phase transition under the bulk 
Gibbs construction. The left panel of FIG. H] shows the mass-central density relations of 
PNS's, and the right panel shows their mass-radius relation. Clearly the high lepton frac- 
tion {Yi = 0.4) enhances the masses because the EOS becomes hard as shown in the left 
panel of FIG. [31 The densities at maximum masses are 3.5-5.0 uq, and the matter of cores is 
in the hadron-quark mixed phase at such densities. The bulk Gibbs construction provides 
a wide density range for the mixed phase as discussed in Sec. II. 
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FIG. 4: Left panel shows the mass-central density relations of PNS's with the hadron-quark phase 
transition under the bulk Gibbs construction, and right panel shows their mass-radius relation. 

Let us move on stable stars with the transition under the Maxwell construction. We find 
that stable stars do not have pure quark phase in their cores just like the case under the bulk 
Gibbs construction, as shown in FIG. O At high lepton fraction (Y/ = 0.4), the maximum 
masses become slightly (3 %) larger than the ones at low lepton fraction {Yi = 0.1). The 
explanation for this is that the phase transition occurs under the density where the effects 
of repulsive nuclear interaction appear on the structures of compact stars. 

There are the central density regions where the gradients of dM/dnB,c are zero, since 
EOS's jump from the hadron phase to the quark phase in these density regions as shown in 
the right panel of FIG. [21 

On the contrary, without the phase transition, the maximum masses at high Yj become 
low clearly (see FIG. [6]). The explanation is that the EOS's without the transition become 
soft at high Yi as discussed in subsection B of Sec. II. 

Maximum masses and the other quantities for each EOS are shown as in TABLE I. The 
label of "bulk Gibbs" ("Maxwell") in the table indicates EOS with the hadron-quark phase 
transition under the bulk Gibbs (the Maxwell) construction, whereas the label of "Shen" 
indicates without the transition, adopting only Shen EOS. 

For the EOS's with the transition, the maximum masses become large for high lepton 
fractions; e.g. for the bulk Gibbs construction, the maximum masses for Yi = 0.1 are in a 
range around 1.95 solar masses. Those for Yi = 0.4 are in another range around 2.10 solar 
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FIG. 5: Same as FIG. HJ but for the Maxwell construction. 
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FIG. 6: Same as FIG. HI but without the hadron-quark phase transition. 

masses which are 8 % larger than the ones for Yi = 0.1. For the Maxwell construction, 
the maximum masses for Yi = 0.4 are 3 % larger than the ones for Yi = 0.1. In this case, 
the typical densities of compact stars are 1 — 2nQ because the maximum densities for stable 
stars are ~ 2nQ. For such densities, the effects of repulsive nuclear interaction at low lepton 
fraction do not appear (see FIG. [2]). 

On the contrary, the maximum masses without the phase transition for Yi = 0.4 are 9% 
smaller than the ones for Yi = 0.1. This is exactly the opposite of the result obtained with 
the phase transition. The following observations help us to understand this difference; the 
maximum densities for stable stars are ~ 4.5nQ. For such densities, the EOS's without the 
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transition become soft at high Yi as discussed in subsection B of Sec. II (see FIG. [2]). 

Let us discuss the evolution of compact stars. We calculate structures of hydro static 
compact stars at a wide range of temperature and lepton fractions. Here, we assume that 
these sets of stars are snapshots of PNS evolutions. Let us change the T and Yi values 
from T = 30 MeV, Yi = 0.4 to T = MeV, Yi = 0.1. The maximum mass difference 
6Mmax between these two cases is ~ O.23M0 (~ O.I3M0) with the transition under the bulk 
Gibbs (the Maxwell) construction. On the contrary, the Shen EOS show opposite behavior 
in maximum masses (i.e. SM^ax ~ — O.IIMq). This difference implies that the EOS's with 
the phase transition become soft after deleptonization although the EOS's become hard 
without the transition. The observations of neutrinos for core collapse supernovae will show 
us which EOS is acceptable in the future. 



B. Unstable stars 



In this subsection, we consider unstable stars which do not satisfy the condition (fT6ll . 

At ub^c/^o ~ 10, all stars are unstable. Such stars collapse to black holes. It is clear 
that a lepton fraction is more important in deciding a stellar mass than temperature (see 
FIG. HUH]). The masses with the transition under the bulk Gibbs (the Maxwell) construction 
are about 1.8 Mq (1.6 Mq) for Yi = 0.1 at any temperature and ns^c/^o = 10, whereas 
they are about 2.0 Mq (1.7 Mq) for Yi = 0.4. This is because the pure quark phase appears 
in these stellar cores. The stiffness of EOS in the NJL model is decided by lepton fraction, 
but not by temperature as described in Sec. II (see FIG. [1]). On the contrary, the masses 
without the transition are about 2.1 Mq for Yi = 0.1, and 2.0 Mq for Yi = 0.4. 

Therefore, these lepton effects on EOS's will be important for some astrophysical phe- 
nomena, such as black hole formations. Our prediction is that trapping of neutrinos will 
suppress the speed of the collapse because the EOS in the NJL model is more stiff. However, 
temperature will not provide such suppression. 



We note that our calculation is only for spherical stars. It is not clear w 
are really stable or not when rotation and magnetic field are also included 



lether the stars 



12 



TABLE I: Maximum masses in FIG. HEl The label of "bulk Gibbs" ("Maxwell") indicates EOS 
with the hadron-quark phase transition under the bulk Gibbs (the Maxwell) construction, whereas 
the label of "Shen" indicates without the transition (adopting only Shen EOS). 



Yi 


T (MeV) 


Mmax/MQ 


R (km) 


nB,c/no 


bulk Gibbs 


0.1 





1.91 


12.5 


4.56 




30 


1.99 


20.7 


3.71 


0.4 





2.05 


12.7 


4.78 




30 


2.14 


19.2 


4.07 


Maxwell 


0.1 





1.94 


13.8 


2.39 




30 


2.02 


24.7 


2.03 


0.4 





1.99 


13.4 


3.38 




30 


2.07 


22.7 


2.80 


Shen 


0.1 





2.31 


12.5 


4.44 




30 


2.38 


16.9 


4.16 


0.4 





2.11 


12.3 


5.00 




30 


2.20 


17.6 


4.46 



IV. CONCLUSION AND DISCUSSION 



We study the structures of PNS's with leptons at finite temperature considering the 
hadron-quark phase transition. To take into account finite size effects, we impose the Gibbs 
or the Maxwell constructions on the phase transition. Our studies show that the pure quark 
phase does not appear in stable stellar cores for all cases. 

The maximum masses with high lepton fraction {Yi = 0.4) become 8 % higher than those 
with low electron fraction [Yi = 0.1) with hadron-quark phase transition under the bulk 
Gibbs construction, because the effects on the NJL model, which stiffen EOS's for high 
lepton fraction, appear in the mixed phase. 

Under the Maxwell construction, the maximum masses with high lepton fraction {Yi = 



13 



0.4) become 3 % higher than those with low electron fraction (YJ = 0.1), because the phase 
transition occurs under the density where the effects of repulsive nuclear interaction appear 
at low lepton fraction. 

On the contrary, without the transition, the maximum masses with high lepton fraction 
become 9 % lower than the ones with low electron fraction, because the hadron EOS which 
we adopt is soft at high lepton fraction. 

Hence, we conclude that the EOS's with the phase transition become soft after delep- 
tonization, although the EOS without the transition becomes hard. 

We also find that lepton fraction is more important to determine structures than temper- 
ature for unstable stars. The explanation for this is that the pure quark phase in the NJL 
model appears in these stelar cores, where the EOS's become hard for high lepton fraction, 
because the chiral restoration for s— quark is suppressed. On the contrary, the EOS's without 
the transition become soft for high lepton fraction as discussed in the last paragraph. These 
behaviors are important for some astrophysical phenomena such as core collapse super novae, 
since the effect of leptons on EOS's changes many dynamical aspects Q, Q, 40, 41, 42|. 

We note that EOS's have many uncertainties. For the hadron phase, we do not take 
into account hyperons in this paper. The effects of hyperons on the phase transition are 
out of the scope of this work and will be presented elsewhere. For quark matter, there are 



uncertainties in the N, 



strong magnetic fie 



PNJL model 



47 



L model itself, such as higher order interaction effects 



43 



45 



4^ and 



461]. Moreover, there are other extended NJL models, such as 



481 . These are open questions for astrophysics and nuclear physics. 
We only use the TOV equation to solve stellar structures. This equation describes stellar 
structures at hydrostatic equilibrium, however it is not clear whether the stars are thermo- 
dynamically stable or not under the same condition. In other words, it does not include 
any effects of thermodynamical convections. Dynamical simulations or linear analysises will 



provide further insight of the effects of those convections. 
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